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We introduce the equilibrium and non-equilibrium statistical hadroniza- 
tion picture of particle production in ultra-relativistic heavy ion collisions. 
We describe the related physical reaction scenarios, and show how these can 
lead to quark pair yield non-equilibrium. Using the SHARE1.2 program 
suite we quantitatively model particle yields and ratios for RHIC-130 run. 
We study how experimental particle ratios can differentiate between model 
scenarios, and discuss in depth the importance of hadronic resonances in 
understanding of hadron production processes. 



1. Introduction 

In strong interaction reaction processes particle production is abundant. 
In the Fermi-Hagedorn statistical model the non-perturbative description of 
particle yields is based on the assumption that particle production is gov- 
erned solely by the size of the accessible phase space. For an introduction 
into the literature and history of the statistical hadronization model we refer 
the reader to Ref . [T] . Further recent developments are discussed in [21 01 • 
We focus in this report on current developments related to RHIC experi- 
mental program, involving the measurement of hadron yields, obtained by 
integrating the produced stable particle spectra. 

* Presented at the XLIV Crakow School of Theoretical Physics 
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The statistical hadronization picture is necessarily an integral part of the 
modeling of this "soft" experimental hadron production data. This is both 
to provide a bottom line to eventual microscopic production mechanisms 
signaling new physics, and to search for abrupt changes associated with a 
phase transformation in which hadrons dissolve into a phase of quarks and 
gluons. 

Hagedorn [1] was first to recognize the importance of hadronic reso- 
nances in statistical particle production. To quantitatively describe hadron 
final state yields, contributions from resonance decays must be taken into 
account. As it turns out, the effect of these resonance decays on particle 
phase space is crucial: As Hagedorn has shown, it can be expected that 
the density density of states of strongly interacting resonances increases 
exponentially with energy, 

p{E) ~ e'^HS (1) 

This has a profound implications for the behavior of matter at high tem- 
perature. It leads to a phase singularity at the temperature ("Hagedorn 
Temperature", Th = ~ 160 MeV) where the canonical partition func- 
tion diverges: 

Z = Y,p{Ei)e-^^/^\T^T^^oo. (2) 

i 

Hagedorn referred to this as the boiling point of hadronic matter. With 
advancement of the quark picture of matter he recognized this singularity as 
being the deconfinement transition, where hadrons cease to be fundamental 
degrees of freedom, and quarks are able to propagate freely throughout the 
system. 

Currently, an intense experimental effort (covering an energy range from 
2 to 200 GeV) is underway to create, and identify, a relatively large volume 
of deconfined quark-gluon matter, the quark-gluon plasma (QGP), and to 
explore its properties. In the following section we describe in turn within the 
statistical hadronization model: particle and quark chemistry in subsection 
12.11 resonance decay chains in subsection 12.21 and we introduce the finite 
resonance widths in subsection 12.31 The reader is referred to Refs. |S] - [20] 
for many further details. In subsection 12.41 we discuss in some detail how 
it is possible that chemical non-equilibrium occurs in a strongly interacting 
system. We than turn to the more phenomenological aspects, and introduce 
the practical aspects of the statistical model in section |31 where we also 
develop strategies for detecting chemical nonequilibrium. In section |1] we 
discus the fitting method and proceed to present some RHIC related results 
obtained with the SHARE suite of programs. 
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2. The statistical particle production model 

2.1. Particle chemistry and yields in the statistical model 

In the grand canonical description the expected baryon and meson yields 
are given by the Fermi-Dirac or, respectively, Bose-Einstein distribution 
functions: 

n{mi,gi;T,Ti) = ni = Qi ( , , (3) 

oo r-p 2 / \ 

= i-2T.i^r-"^^^Kj^). (4) 

n=l 

In the upper signs refer to fermions and the lower signs to bosons, re- 
spectively. The second form, Eq. (j^, expresses the momentum integrals 
in terms of the modified Bessel function This form is practical in the 
numerical calculations. The series expansion (sum over n) converges when 
TjC"™''/"^ < 1, below the Bose-Einstein condensation limit. Consideration 
of this condition is required only for the pion yield in the range of parame- 
ters of interest. Here, the index 'i' labels different particle species, including 
hadrons which are stable under strong interactions (such as pions, kaons, 
nucleons or hyperons) and hadron which are unstable {p mesons, A (1232), 
etc.). rrii is the particle mass, while the quantity gi = (2Jj -|- 1) is the spin 
degeneracy factor; we will distinguish all isospin sub-states separately. 

Tj is the fugacity factor. Since particle anti-particle pairs (ii) can be pro- 
duced by particles without conserved 'charges' (baryon number, strangeness, 
etc.), 

i + i ^ n X 7, (7, TT, ...(n — objects with no conserved charges) (5) 

a scenario where statistical production arises as a limit of kinetic evolution 
requires that the chemical potentials of particles and antiparticles associated 
with the conserved charge be opposite. Thus the fugacity Tj comprises the 
well known factor A^^ = e^^/'^ , where is the chemical potential associated 
with a conserved quantum number. 

For a dynamically evolving system yield of particles cannot in general 
be in chemical equilibrium. This is accommodated by introducing a phase 
space occupancy factor 7: 

T^i=7iAj, Tj = 7jA-^ (6) 

7j = 1 corresponds to chemical equilibrium. 7^ < 1 means that quantum 
number j is under-saturated, usually because time is needed to build it up. 
An isolated system kept long enough "in a box" would generate more of 
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jj particle-antiparticle pairs at the expense of the internal energy. If this 
equilibrated box undergoes isentropic and rapid expansion {VT^ = Const, 
for massless particles) , 7^- > 1 wih generahy result as the particle abundances 
do not have time to re-equilibrate Given time in the "large box", pairs 
would annihilate and the system would again relax to an equilibrium state. 
In subsection 12.41 we will explore these scenarios more thoroughly, and in 
section |31 we will argue that equilibrium assumptions can and should be 
tested against experimental data. 

Each hadron can be a carrier of several conserved charges contained in 
its valance quarks. The most general fugacity Tj for hadron j is given by a 
suitable product of A and 7 factors: 

= A,. (A,7.)'^' (A.7.)'^^' iKlcf' (A.-7.)^^' (A.-7c-)^^' , (7) 



'3 



where 



and 



Ag — Aq^, As — Ag^, Ac — Ag^, (8) 



Iq = 7q, Is = 7s, 7c = 7c- (9) 



Here, A^^, and A'^^ are the numbers of light {u, d), strange (s) and charm 

(c) valance quarks in the j'th hadron, and iV|, and are the numbers 
of the corresponding valance antiquarks in the same hadron. 

The quark and hadron chemical potentials as used in e.g. Refs. [T^ 1^ 
I23j are related by simple algebraic conditions arising from the valance quark 
content and definitions of the baryon and hyperon charges of particles given 
in the historical context. Considering the case of chemical equilibrium limit 
7q = 7s = 1 and in absence of charm (A*"* = = 0): 

rr = exp ( EltE±^ltl±Iil^] , (10) 



where Bi, Si, and are the baryon number, strangeness, and the third 
component of the isospin of the ith particle, and ^u's are the correspond- 
ing chemical potentials. In this case, the two formulations are related by 
equations: 



A^^? = A,. = e^3>^3/^. 



(13) 
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We note that strangeness and hyperon number S" of a particle are by his- 
torical convention opposite in sign to each other, 

Since we use the symbol S also for entropy, henceforth we will refrain from 
mentioning the hyperon number. 

2.2. Resonance decays 

In the first instance, we consider hadronic resonances as if they were par- 
ticles with a given well defined mass, e.g., their decay width is insignificant. 
All hadronic resonances decay rapidly after freeze-out, feeding the stable 
particle abundances. Moreover, heavy resonances may decay in cascades, 
where all decays proceed sequentially from the heaviest to lightest particles. 
As a consequence, the light particles obtain contributions from the heavier 
particles, which have the form 

ni = b2^i ...bN^N-inN, (14) 

where b^^k-i combines the branching ratio for the k ^ k — 1 decay (ap- 
pearing in j24j ) with the appropriate Clebsch-Gordan coefficient. The latter 
accounts for the isospin symmetry in strong decays and allows us to treat 
separately different charged states of isospin multiplets of particles such as 
nucleons. Deltas, pions, kaons, etc. For example, different isospin multiplet 
member states of A decay according to the following pattern: 



A++^7r+-hp, (15) 

A+^i(7r++n) + ^(7r0+p), (16) 

AO-i(vr-+p) + ^(vrO + n), (17) 

A^^vr^+n. (18) 



Here, the branching ratio is 1 but the Clebsch-Gordan coefficients introduce 
another factor leading to the effective branching ratios of 1/3 or 2/3, where 
appropriate. 

A decay such as A — > vrA^ is easy to deal with, since only one de- 
cay channel is present and the branching ratio is well known. However, 
for most of the resonances several decay channels appear and the partial 
widths (product of branching ratio with total width) are in general less well 
known, usually classified as dominant, large, seen, or possibly seen, |24jl. 
Any statistical hadronization calculation where high-lying massive states 
are relevant should carefully consider the dependence of this uncertainty in 
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its results, especially considering the increasing resonance degeneracy with 
mass (221 • 

Weak hadron decays complicate the physics of hadron yield and spec- 
tra study considerably. The weak interaction life span for strange hadrons 
(excluding K^) ranges in 2-10 cm/c, which values are subject to the time 
dilation effect in the decaying particle rest frame. In general these weak de- 
cay feed-down (WDFD) occurs both near and far enough from the primary 
interaction vertex so that a little known fraction of decays is counted, or not 
counted, depending on the type of the experiment. For experiments mea- 
suring protons at RHIC, the WDFD by hyperons is large, in the vicinity of 
50-60%. For antiprotons at AGS energies the contamination by antihyperon 
decay is dominant. Similarly, the yield of pions includes a large fraction of 
Kg decay products and some hyperon decays. Sometimes experimental re- 
sults include corrections obtained in simulations of WDFD, stripping the 
reported yields of the WDFD contamination. These corrections are specific 
to each experiment and even to the experimental cuts made in the particu- 
lar analysis. Such corrections introduce additional systematic uncertainties 
since model dependent assumptions about the primary particle spectra must 
be made. Often enough given these complications, the reported results are 
not corrected for WDFD contamination, leaving us with an interesting mea- 
surement which is hard to compare to a theoretical model. 



For a particle i which has a finite and significant (thus hadronic) decay 
width Fj, the thermal yield is more appropriately obtained by weighting 
Eq. ((31) over a range of masses in order to take the mass spread into account: 



The exponential thermal weight n(M, g^; T, Tj) is asymmetric around the 
value of the mass of the particle M = rrii and when Fj and T are of compa- 
rable magnitude we find that this in general increases the yield as compared 
to the Fj — > limit. But that is not the end of the story. 

The Breit-Wigner distribution if used with energy independent width 
means that there would be a finite probability that the resonance can be 
formed at unrealistically small masses. Since the weight involves a thermal 
distribution n{M, gi;T,Ti) which would contribute in this unphysical do- 
main, one has to use, in Eq. (fTH)). a more physical, energy dependent width. 
The dominant energy dependence of the width is due to the decay threshold 
energy phase space factor, dependent on the angular momentum present in 



2.3. Finite resonance width 




for Fi ^ 0. (19) 
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the decay. The exphcit form can be seen in the corresponding reverse pro- 
duction cross sections [26n2Zj- The energy dependent partial width in the 
channel i ^ j is to a good approximation: 



M 



'»J ^ 2 



for M > m, 



(20) 



Here, rriij is the threshold of the decay reaction with branching ratio bi^j. 
For example for the decay of i := A"*"^ into j := p -\- vr^, we have niij = 
TUp + , while the branching ratio is unity and the angular momentum 
released in decay is lij = 1. From these partial widths the total energy 
dependent width arises, 



F, ^ri(M) = ^ri_,(M). 



(21) 



For a resonance with width, we thus have replacing Eq. (|19() : 



1 r°° 
'^i=wT. dMn{M,gi;T,T,)- 



(M) 



(M-m02 + [r,(M)]2/4' 
and the factor (replacing 2tt) ensures the normalization: 



E 



dM 



(M-mi)2 + [r,(M)]2/4- 



(22) 



(23) 



The resonance yields with widths in general vary considerably when com- 
pared to the yields without widths. One finds now both enhanced and 
suppressed yields as compared to the limit Fj — 0. 



2.4- Chemical nonequilibrium 

Consider a large "pot of boiling quark-gluon soup" . In the statistical 
hadronization approach hadrons 'evaporate' with an abundance correspond- 
ing to the accessible phase space. However, the quark equilibrium in the pot 
in general precludes that that just after formation, the evaporated hadrons 
are in chemical yield equilibrium. In particular, it has long been understood 
j28j that the quark phase space and the hadron phase space are different. 

It is therefore natural to expect that, if chemical freeze-out is near a 
phase transition in which hadrons are formed by quarks recombining ac- 
cording to hadronic phase space, the hadron gas will be in thermal but not 
chemical equilibrium. Entropy conservation suggests that, if the transition 
is rapid (so the system can not accommodate e.g. a smaller number of 
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available degrees of freedom through expansion), particle phase space will 
be oversaturated significantly above the point at which detailed balance is 
achieved [5^]. We could therefore expect ^q^s > with 7^ approaching the 
point of B-E condensation (e™'"/^^). 

Another physics mechanism could lead to the inapplicability of chemical 
equilibrium to strange quark abundance. Strangeness is produced relatively 
slowly through kinetic reactions (between quarks or hadrons). It is natural 
to expect the abundance of strange and antistrange particles to approach 
equilibrium asymptotically from below. If the system's lifetime is too short, 
strange particle abundance is expected below chemical equilibrium, 7s < 1 

IHOI. 

Strangeness suppression in a small system {e.g. pp,pA), as compared 
to grand canonical yield Eq. can also be understood within statistical 
model as due to phase space reduction in the canonical statistical model |^ . 
When strangeness enhancement in AA collision is reported based on com- 
parison with respect to these smaller systems, an artificial enhancement is 
expected |32| . This artificial enhancement effect depends in a specific fash- 
ion on volume, particle strangeness, and reaction energy. This mechanism 
does not describe the AA data in detail Recent experimental results 
confirm this. Both the predicted rapid approach to the grand canonical limit 
as function of the number of participants, and the enhancement increase of 
e.g. the H, S with decreasing energy disagree decisively with experimental 
40 GeV and 158 GeV SPS results j^i^SilSl!- With these experimental 
results the effect of strangeness enhancement in AA reactions is confirmed 
to be the result of kinetic production processes. 



3. Statistical hadronization model 

3.1. Model- data comparisons 

In the most general version of the statistical model, (excluding in this 
consideration charmed quarks), we have six independent parameters which 
have to be determined by experimental data: T, A^^s, 7^^^, A/3, and an over- 
all volume parameter normalizing particle yields. This volume V is pro- 
portional to the magnitude of the transverse geometric area of the reacting 
system determined by the size of the nuclei and the event trigger, and it 
depends on the model dependent longitudinal rapidity acceptance relation 
with the longitudinal volume extend. This dependence of V on non-trivial 
and also experiment-specific factors makes it tempting to eliminate it in a fit 
of particle ratios. For yield studies which include the full particle multiplic- 
ity, the absolute yields can, however, be fitted and the volume parameter is 
determined by the trigger condition which fixes the number of participants. 



zakopanevl.3 printed on February 8, 2008 



9 



One of the remaining 5 parameters, /a, is fixed by the proton to neutron 
asymmetry in the reaction. No matter how uncertain we are about the 
longitudinal dynamics, the ratio ((d) — {d))/{{u) — (u)) (where (d) refers to 
the number of valance d-quarks, etc.) is fixed and we accommodate this by 
fixing the charge to baryon ratio of the colliding nuclei. Another constraint 
is strangeness conservation, that is {s) = (s). When implementing this 
condition one generally fixes the value of Ag. Here there are, however, some 
caveat. 

First, since strangeness decays weakly, the final strangeness yield is, 
at finite baryon number, asymmetric. Considering that there are more 
strange baryon weak decays than strange antibaryon weak decays, the true 
weak decay correction situation is that the final state has (s) < (s). In 
consequence of experimental corrections which tend to counterbalance this 
effect one should perhaps not impose exact strangeness conservation but 
allow a measure of uncertainty when implementing strangeness conservation. 
Secondly, because of baryon asymmetry, the rapidity distribution of strange 
and antistrange baryons needs not to be exactly identical. Given that the 
baryon number at RHIC is expected to be mostly in fragmentation regions, 
this is where one could expect some excess of (s) > (s). 

We obtain the fit parameters from the experimental by minimizing x^- 
When fi{T, A, 7) is the model prediction for particle i, ni is the experimental 
value, and Astat,sys are the statistical and systematic errors. 



The best fit parameters T, Aj, 7^ are those which minimize while respect- 
ing constraints described above. 

The previous section makes it clear that any fit of statistical model pa- 
rameters to data should, in principle, allow for the possibility that some, 
or all quantum numbers are not in chemical equilibrium. In other words, 
equilibrium should be an eventual result of the study, rather than it's as- 
sumption. However, fixing the non equilibrium parameters to a tacit value 
such as jq^s = 1 reduces seemingly (that is by assumption) the number of 
parameters, which increases n, the number of degrees of freedom (d.o.f) 
in the fit. To compare fits which use jg^g with those which fix it in case 
that n is not very large, it is necessary to consider a comparison standard 
independent of the value of n. Such a standard is provided by statistical 
significance P shown in Fig. ^ Ref. [21] . Each line is drawn for a fixed 
value of P [%] defined as the likelihood, given the validity of the model, 
and a purely Gaussian error source, of finding the corresponding y^^ jn^ or a 
smaller value. 




(24) 
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10 20 30 40 50 

Degrees of freedom n 

Fig. 1. Lines of constant statistical significance P [%] as a function of the number 
of degrees of freedom n and ^ jn |24 |. 



As can be seen P reaches an asymptotic value for n — > oo, where it 
rapidly approaches 100% for ^ jn < 1 and rapidly approaches otherwise. 
However, the typical number of degrees of freedom in statistical model fits is 
considerably below the asymptotic limit, making the statistical significance 
dependent on n. For this reason perfectly acceptably-looking graphs can 
have an unacceptably small statistical significance. For example, for a fit 
with e.g. 5 measurements and 3 parameters to have P = 90% it must come 
up with ^ 0.1, that is the fit must hit the center of each measurement 
to be considered based on a good physical model. The range X'^/^ ~ 1 
corresponds in this case to P = 30%, which means that most likely the 
physical model used is false or the data inaccurate or too precisely stated. 

Statistical significance, per se, is not a "proof of validity" of a model. 
It is, however, a discriminating force. If a fit has a 'small' statistical signif- 
icance, well below 60%, the fit's validity and/or experimental data should 
be questioned, e.g. are there unaccounted for systematic errors? Are par- 
ticular data points spoiling the fit, and if so, why? Thus study of P must 
supersede all "looks nice" arguments. And if a markedly higher statistical 
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Table 1. RHIC-130 GeV hadronization parameters 



parameter 


r = 


vary 

r [211 


79.S 

r = 


= 1 

r^o 


Iq = 1, 7. 

r = 


. varies 
r ^ 


T [MeVJ 


133 ± 10 


135 ± 12 


158 ± 13 


0.157 ± 15 


152 ± 16 


153 ± 23 


lO^A, - 1) 
102(A«i-l) 


708 ± 342 


703 ± 337 


735 ± 390 


730 ± 382 


724 ± 373 


721 ± 363 


3.132973 


3.203555 


2.636207 


2.788897 


2.95346 


3.00848 


7? 


1.66 ± 0.013 


1.65 ± 0.030 


1 


1 


1 


1 


7s 


2.41 ± 0.61 


2.28 ± 0.46 


1 


1 


1.17 ± 0.30 


1.10 ± 0.25 


104(Az3 - 1) 


30 ± 305 


28 ± 293 


59 ± 564 


53 ± 508 


64 ± 525 


59 ± 481 




fit relevance 


N — p = n 


16-5 


16-5 


16-3 


16-3 


16-4 


16-4 




0.4243 


0.4554 


1.0255 


0.8832 


0.6067 


0.7301 


P[%] 


94.61 


93.07 


42.25 


57.05 


83.85 


72.32 



significance can be acfiieved by varying from a tacit value a physically moti- 
vated fit parameter, this can be taken as evidence that the physical scenario 
underlying this parameter is relevant. 

3.2. Sample of hadron yield fit results 

The table ^ and Fig. ^ were also presented in [HH] ■ Similar results are 
found in fits to SPS data [HHI- We note that Tabled shows that strangeness 
under-saturation models appear to be ruled out at RHIC-130. The fitted 
parameter for a non equilibrium fit is well above unity, and it is slightly 
above above unity in fits where 7^ = 1. This latter observation is corrob- 
orated by another fit of RHIC data We recall that SPS and AGS 
results are best fitted with ^sMq < 1 HH- We see in tabled that chemical 
equilibrium model fits yield relatively low statistical significance P, which 
is increasing as we step up the non equilibrium, introducing first 7^ 7^ 1 and 
than 7q 7^ 1 . 

Even though there are some serious differences in confidence level, the 
three models look nearly indistinguishable in Fig. ^ Fig. ^ shows that at 
~ 10 degrees of freedom, one needs a x^/n DoF considerably below 1, which 
requires that the average data point appears within the error bar. 

Given the present data quality, it is not possible to determine with pre- 
cision values of 7^^^ , and hence to differentiate full non-equilibrium from 
strangeness non-equilibrium. We can try constraining the fit parameters 
further by examining sensitivity of data point to them in detail. Figures |31 
|1] and IS] aim to to that. To prepare them, we have performed the best fit 
at a given fixed T for Fig. |31 for |1] and 7^ for |SJ varying all the other 
parameters. The model data points are than graphed against the parame- 
ter in question to show sensitivity to the parameter in this observable. It is 
immediately apparent that, in models with chemical non-equilibrium, pa- 
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n /n z^/Z 
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K77i" n7n- 



TVA k7K" 

K7/ 



*o :i 
K /K" 



A_ 

(|)/h" 



Fig. 2. A comparison of 130 GeV Au-Au data and the three different statistical 
models - all results 'look good', even if statistical significance varies considerably, 
see bottom of tabled 



rameter correlation makes it possible to model some particle ratios over a 
wide range of any parameter. For instance, Fig. |31 shows that once temper- 
ature reaches ~ 130 MeV, it is possible to model ratios such as K/tt and 
'E/Q using any temperature within the range. It is only when several similar 
ratios are fitted together that parameters become disentangled. 

Figs EES also show that some data points are better at disentangling fit 
parameters than others. Fig. El shows that resonances (here K*) are ideal 
for constraining the freeze-out temperature. If compared to a stable parti- 
cle whose quark composition is identical, all chemical factors (A, 7) cancel 
out independently of equilibrium, so the ratio is a function of temperature 
and flavor content alone. When the K*/K ratio has been measured to a 
greater precision, it will be possible to definitely constrain the freeze-out 
temperature. This could be sufficient to distinguish an equilibrium scenario 
from one where 7^ is needed to explain hadronic abundance. 

Constraining chemical potentials, in particular to ascertain directly if 
'jq^s 7^ 1 are present, is more difficult, since these parameters are strongly 
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' • ' ^ ' • ' — ^ • ' • ' 

0.1 0.12 0.14 0.16 0.18 0.2 

T [GeV] 



Fig. 3. Sensitivity of a measurement of K*'^/K~ , K~ /7t~ , Vl/'E.~ on Temperature 
T, other statistical hadronization model parameters have been optimized to best 
fit the other experimental RHIC-130 data. 



correlated with temperature. Thus measurement of resonances also fixes 
this issue, in an indirect way. In addition, as Fig. El shows, 7^ can be 
constrained by measuring ratios of particles of similar mass, but different 
strangeness content, such as H and A, while 7^ (Fig. 0]) is sensitive to ratios 
between baryons and mesons. At the moment, the error bars for these data 
points are too large to allow disentangling T, 7^ and 7<j. 

4. How to SHARE 

In the following years, the high-statistics RHIC date, beginning with 
200 GeV run will be thoroughly analyzed, leading to better determination 
of the yields of many stable particles and resonances. As shown in the 
previous section, these measurements should be sufficient to distinguish be- 
tween models described above. In addition, SPS and eventually the future 
GSI facility, as well as higher energy LHC measurements promise an exten- 
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Fig. 4. Sensitivity of a measurement of /p, /p^ , /tt^ on the value of 7^, 
the chemical light quark pair occupancy parameter, other statistical hadronization 
model parameters have been optimized to best fit the other experimental RHIC-130 
data. 



sive energy scan, allowing a full assessment of the applicability of different 
statistical models in a wide range of ^/s, as well as the detection of any 
systematics in the fit parameters. 

We have argued that, if the purpose of statistical analysis is to observe 
an eventual phase transition, any assumptions about chemical equilibration 
have to be tested against experimental data. We have also argued that 
statistical hadronization calculations require a very detailed input of the 
hadronic spectrum, in particular in regard to the uncertainty in degeneracy 
and decay patter of the higher-lying states [23 as well as experiment- 
specific weak-decay feed-downs. 

To get the correct physics out the experimental data, and to ascertain 
the physical significance behind the statistical model, these approaches need 
to be compared with experimental data using a standard procedure. We 
have proposed an open source general model for: Statistical HAdronization 
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Fig. 5. Sensitivity of a measurement of Q/'ET , S^/A, 'Br /<j) on the value of 7^, the 
chemical strange quark pair occupancy parameter, other statistical hadronization 
model parameters have been optimized to best fit the other experimental RHIC-130 
data. 



with REsonances (SHARE) This statistical hadronization package 

is suitable for comparative analysis of heavy ion collisions in a variety of 
systems, energy ranges, and equilibrium assumptions. A cross check with 
SHARE requires that: 

1. the resonance decay tree be the same for the different models under 
consideration, so as all systematic effects in the fit parameters and the 
quality of the fit deriving from choice of resonance decay are under 
control, 

2. experiment-specific weak decay acceptances should be treated sepa- 
rately for each experimental data-point, to minimize systematic errors 
on fit parameters deriving from experimental effects, 

3. each model's ability to fit data should be compared in a way that takes 
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into account the different number of degrees of freedom within each 
model. Statistical significance is a comparison criterion which satisfies 
this requirement, 

4. the fitting program should be able to fit each model's fit parameter, 
or constrain it using conservation laws (if applicable), 

5. the fitting program should be able to evaluate the fit's sensitivity to 
each parameter checking for false minima, etc. and employing 
profiles and contours, 

6. the fitting program should be able to test the fit's sensitivity to each 
data point. This is invaluable both to determine the consistency of 
the fit, and to detect systematic violations indicating possible novel 
microscopic formation mechanisms, 

7. to better assess the physical soundness of the statistical picture, bulk 
quantities (total number of negative/positive/neutral particles, en- 
tropy, pressure, energy density etc. should also be calculated. 

All these features are part of the package SHARE |121- Moreover we of- 
fer three different evaluation of particle multiplicities (SHARE FORTRAN, 
Mathematica SHARE, and web based SHARE - derived from the FOR- 
TRAN version of SHARE). SHARE calculates yields, ratios and bulk quan- 
tities in terms of thermodynamic parameters, Fortran version performs fits, 
and explores parameter sensitivity to data (x^ profiles, contours etc.). The 
resonance decay trees, and weak corrections are allowed form and Fortran 
version requires in the fit the experimental ratios as input. 

The reader interested in this field is encouraged to use SHARE in study 
by means of statistical hadronization models new experimental data. The 
program is available online, at 

http : //www . physics . arizona . edu/ ~torrieri/SHARE/ share . html 
http : //www . if j . edu . pl/Dept4/share . html 

We hope that this effort will contribute in achieving a better understand- 
ing of the applicability and physical significance of the statistical particle 
production picture. 
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